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A motile cell, when stimulated, shows a dramatic increase in the activity of it's 
membrane, manifested by the appearance of dynamic membrane structures such as 
lamellipodia, filopodia and membrane ruffles. The external stimulus turns on mem- 
brane bound activators, like Cdc42 and PIP2, which cause increased branching and 
polymerization of the actin cytoskeleton in their vicinity leading to a local protru- 
sive force on the membrane. The emergence of the complex membrane structures 
is a result of the coupling between the dynamics of the membrane, the activators 
and the protrusive forces. We present a simple model that treats the dynamics of a 
membrane under the action of actin polymerization forces that depend on the local 
density of freely diffusing activators on the membrane. We show that, depending 
on the spontaneous membrane curvature associated with the activators, the result- 
ing membrane motion can be wave-like, corresponding to membrane ruffling and 
actin-waves, or unstable, indicating the tendency of filopodia to form. Our model 
also quantitatively explains a variety of related experimental observations and makes 
several testable predictions. 



I. INTRODUCTION 

Various types of directed cell motility are driven by the polymerization of an actin network 
that exerts a force on the cell membrane, pushing it forward. During cell motility, the leading 
edge of the cell exhibits a range of dynamic structures such as lamellipodia, filopodia and 
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membrane ruffles (Small, 2002;Pollard, 2003). These dynamic surface patterns of moving cells 
are usually observed to have lengthscales in the 1/im range, and appear in many different 
cell types (Gungabissoon,2003). The lamellipodium is a flat, disk-like extensional structure 
generally occurring at the periphery of spreading cells while the filopodia or micro-spikes 
are actin rich needle-like structures seen generally as extensions of the lamellipodium. To 
generate movement, the cells use precursor contacts found in membrane ruffles, or on the 
underside of filopodia, which can help form adhesive contacts. The lamellipodia are generally 
extruded in the direction of a strong signal such as a chemoattractant that induces cell 
migration. The extracellular stimulus turns on certain membrane bound activators that in 
turn activate a series of proteins that trigger actin branching and polymerization leading 
to a directed and regulated protrusive force. The interplay between the dynamics of the 
activators, the protrusive forces generated by the actin polymerization and the membrane 
dynamics results in the rich variety of dynamic structures described above. This kind of 
actin based motion is ubiquitous with examples ranging from the chemotaxis of macrophages 
to the movement of metastatic tumor cells. Since cell motility depends so crucially on the 
formation of these dynamic membrane structures, it is imperative to understand the origin 
and dynamics of these structures. 

There have been a number of theoretical approaches to the problem of cell motility driven 
by actin polymerization (Mogilner,1996;Bottino,2002;Grimm,2003;Stephanou,2004). These 
studies focus on the inter-relationship between the dynamics of actin polymerization and 
the protrusive forces generated which lead to propulsion. However, one important aspect 
that has been neglected so far is the crucial role played by membrane bound activators 
(henceforth referred to simply as activators or membrane proteins), and in particular the 
thermal density fluctuations and the spontaneous membrane curvature associated with the 
activators. Activators that have been well studied include the Rho family GTPase Cdc42 
and the membrane phospholipid PIP2 (Higgs,2001). These activators bind to, and activate, 
WASp/Scar family proteins by inducing a conformational change. The WASp/Scar pro- 
teins in turn activate Arp2/3 which is directly responsible for generating new branches in 
the actin network (Higgs,2001;Blanchoin,2000;Carlier,2003). Fluctuations in the density of 
the membrane activators can thus directly lead to fluctuations in the actin branch density 
and hence the protrusive force. Also important is the fact that both Cdc42 and PIP2 have 
been shown to induce and sense membrane curvature by binding to BAR domain proteins 
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and epsin respectively (Habermann,2004;Ford,2002). In this paper we propose an approach 
that takes into account density fluctuations and membrane curvature associated with the 
membrane activators. The dynamics of the actin polymerization, and its dependence on 
the relative concentration of actin and supporting proteins, has been calculated by Carlsson 
(Carlsson,2003). This gives us the steady-state velocity of the advancing membrane and 
actin gel, as a function of the branching rate which is directly proportional to the mem- 
brane density of the activators. We therefore separate the in-cell dynamics of the actin 
polymerization, from the in-membrane dynamics of the activators. The diffusion and spon- 
taneous curvature associated with these membrane activators will determine the time and 
length-scales of the dynamic patterns. This allows us to write a simple model and arrive 
at analytical expressions, while still preserving the rich variety of dynamical behavior that 
is observed. The merits of such an approach, apart from the knowledge gained concerning 
dynamical models, lie in its quantitative and testable predictions for cell motility in vivo. 

Our work draws on previous models of instabilities in active membranes 
(Chen,2004;;Sankararaman,2004;Manneville,2001;Ramaswamy,2000;Prost,1998), some of 
which we have combined here into a simpler form. The previous analyses are different 
in essential ways from the present study. They consider the case where the active membrane 
proteins are ion-pumps and therefore solve for the fluid circulations, which is not neces- 
sary here. The first analysis (Chen, 2004) considers a system close to a phase separation 
transition, which is not a constraint in our work. The possibility of wave-like propaga- 
tion for negative spontaneous curvature of the membrane proteins, was also not considered 
(Chen, 2004;Sankararaman, 2004). The second analysis (Manneville,2001;Ramaswamy,2000) 
does find propagating modes, but does not specifically relate them to the spontaneous cur- 
vature of the membrane activators. Our model is therefore an application of these previous 
studies to the case of actin-driven cellular motility. It allows us to describe, in a very detailed 
and transparent way, the physics of this system in terms of an active membrane. 

II. MODEL 

We now introduce our model for the dynamics of the coupled system consisting of the 
membrane, the activators and the actin polymerization induced forces. The problem we 
wish to solve is shown schematically in Fig.l. We have an average areal density of activating 
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proteins, n , which induces actin polymerization such that the membrane moves forward at a 
velocity given by: v = An ~ 1— O.lum/sec, where A is a coefficient that depends on various 
factors, such as concentration of actin monomers, temperature etc., and no ~ 2 • 10 15 m -2 . 
The time evolution of the density of these activating proteins is described by a diffusion 
equat ion ( Manneville ,2001; Ramaswamy, 2000) 

— = DV 2 n - AnHW^h + V • f n (1) 
at 

where D ~ lum 2 /sec is the in-membrane diffusion coefficient of the proteins, k is the bending 
modulus of the membrane (typically k ~ 10k B T), A = D/x is mobility of the proteins in the 
membrane where x — k B T is the effective in-plane compression energy of the proteins, H is 
the spontaneous curvature of the membrane proteins, and the thermal noise force satisfies 
the following correlation: (f n (r,t)f n (r',t')) = 2n Ak B T5(r — r')5(t — t'). The first term 
describes the free diffusion of the activators on the membrane, while the last term captures 
the effect of thermal noise. Evidence for free diffusion of the membrane proteins that activate 
the polymerization of actin, appears in (Gerisch,2004). The second term takes into account 
the coupling between the spontaneous curvature of the activators and the local curvature of 
the membrane. Here h refers to the coordinate that measures the normal displacement of 
the membrane from a flat reference plane. 

The membrane deviation from flatness obeys the following equation of motion 

^ = - J dr'A(r - r')KV%{r') + An (2) 

for a free, flat membrane. Here the first term is simply the response of a membrane that is 
surrounded by a fluid, characterized by the hydrodynamic interaction kernel A. After Fourier 
transforming into g-space, the hydrodynamic interaction kernel is given by A(q) = l/Arjq, 
and the response of the free membrane is u q = Kq 3 /4ri, where 77 is the viscosity of the 
surrounding fluid. We note that if the edge of the membrane is highly curved, then the 
membrane response is different. We take the response in this limit to be: u qt i ~ Kq/ir/d 2 , 
where d is the local radius of curvature at the membrane edge (Fig.l), assuming that the 
hydrodynamic interaction (Oseen kernel) remains the same as that for the flat membrane 
(Zilman,2002). A similar form for the response of a flat membrane occurs when tension is 
dominant and is given by u q;t = aq/2i], where a is the effective surface tension. This regime 
occurs for wavevectors: q < q t = Jcr/n. In the analysis below, the results using either u Qj i 
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or u q> t are interchangeable, by making the transformation: k/2cL 2 <-> a. The second term 
describes the action of the actin polymerization induced forces, whose effect we model by the 
addition of a velocity to the membrane. This velocity is taken to be directly proportional to 
the local density of activators, with a proportionality constant, A, as described above. The 
two equations (eqs. II 1211 form a coupled set that completely describes the dynamics of our 
system. We now describe, in detail, the assumptions that we make in our model. 

A. Discussion of assumptions 

In this paper, we assume that the membrane proteins do not bind to the actin network, 
and that the diffusion coefficient (in eq^) is homogeneous. The diffusion coefficient, D, 
may be treated as an effective value, which takes into account the average (uniform) effect 
of actin interactions with the membrane proteins. This is an approximation, since the 
diffusion coefficient is likely to decrease when the density, n, increases, due to a "crowding 
effect" (Almeida, 1995). It may also depend on the actin density. Future work may include 
the dependence of the diffusion coefficient on the actin density, membrane curvature and 
the dynamics of protein-actin binding/unbinding. The binding of the membrane proteins 
to the actin network that they nucleate, introduces an effective attraction between them, 
which may drive phase separation. An effective protein-protein interaction will introduce a 
term of the form: J(r — r')n(r)n(r') , to Eq^ This possibility is deferred to future studies. 

We use the term "activator-proteins" in the most general sense: it stands for membrane 
proteins that trigger actin polymerization, branching or bundling, that produces in turn a 
protrusive force on the membrane. In addition to the activators mentioned above, there 
are, for example, the VASP membrane proteins that recruit fascin, which cross-links actin 
filaments into bundles. These bundles can push more effectively on the membrane, and 
produce filopodia (Biyasheva,2004). We consider a single species of membrane activator, 
which is also constantly in its "on" state. A more detailed description could allow for the 
kinetics of the turning-on of these activators, which is a process influenced by the presence 
of other proteins, such as external chemotactic signals or binding of other cellular proteins, 
the average density n and membrane curvature. Nevertheless, in this paper, we wish to 
investigate the dynamics that arise from the simplest model first. 

In EqEJ we neglect the thermal fluctuations of the membrane, which are usually much 
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smaller than the motion due to the actin polymerization (see Section IV). The thermal 
fluctuations of the membrane and the actin-induced motion described above are incoherent 
(decoupled) , so that they simply add to the overall mean-square height fluctuations. We 
further expect the thermal fluctuations of the membrane to be almost eliminated when the 
membrane is being pushed by the actin network, since any membrane motion that is not 
synchronized with the actin polymerization, such as the thermal motion, will be reduced 
to negligible values due to the large bulk modulus Y of the actin gel, giving mean-square 
height fluctuations: (h 2 ) oc ksT/Y. The membrane response, u q , which we used above, for 
a membrane that is pushed by the actin polymerization, describes the dynamics of the fluid 
flow outside the cell, and does not need to include the elasticity of the actin network. For the 
same reason, membrane motion that arises directly from the spontaneous curvature of the 
membrane proteins (Chen, 2004;Manneville,2001;Ramaswamy,2000;Sankararaman, 2004), is 
also negligible since it too is incoherent with respect to the motion due to actin polymeriza- 
tion. 

In EqElwe also assume that there is a linear relation between the density of activators, 
n, (and therefore of activated Arp2/3 protein) and the forward velocity of the membrane. 
That this assumption is valid for low densities and velocities (i.e. velocities low with respect 
to the saturation polymerization velocity: v p ~ 1/im/sec) has been shown explicitly within 
the context of a model that considers an obstacle driven forward by a polymerizing actin 
network with a spatially homogeneous branching rate (Carlsson,2003). Above a very small 
cut-off branching rate, the forward velocity increases linearly with the branching rate and 
saturates to a maximum at high branching rates. This behavior is generic and not crucially 
model dependent. One can also reach the same conclusion from a continuum perspective. 
At low densities the branching of the actin gel increases linearly with the activator density 
(Blanchoin,2000). This means that the bulk modulus of the gel will aslos be a linear function 
of the activator density: Y oc n. The relation between the velocity normal to the plane of 
the membrane and the modulus of the pushing gel is given by (Gerbal,2000) 

- = (3) 
v p 1 + 

At low gel densities, Eq|H] implies: v/v p ~ Y(Sb/F ext ), where F ext is the external drag force 
and Sb is the local area of the membrane which is pushed by the actin gel. Thus we do 
find a linear relation between the velocity and the density n. At high densities the velocity 
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saturates at v p , and our analysis breaks down. 

Note that our description of an imposed velocity, v(n) = An in Eq|21 due to the 
actin polymerization, is different from that of an imposed force condition F ac u n (Man- 
neville,2001;Ramaswamy,2000;Gov,2004). The imposed velocity condition is natural if it is 
determined by the dynamics of the actin polymerization in the lamellipodia. This condition 
also applies if there is a roughly constant drag force due to the actin gel itself (Gerbal,2000). 
On the other hand, if the motion is determined by the action of the drag force of the sur- 
rounding fluid, then it is more natural to keep the imposed force condition, and substitute: 
A -> F actin n/r]q, in Eq0 

Another assumption implicit in EqJ21 is that changes in the membrane density of the 
activators translates instantaneously into changes in the force with which the actin gel is 
pushing the membrane. This is not strictly true, and we now wish estimate the time lag for 
this process. First there is the chemical time for Arp2/3 activation. This is of order 1msec, 
which translates to membrane density fluctuations of length-scale lOnm, which, in turn, is 
much shorter than the typical mesh size of the actin gel (~ 50nm). So we may neglect this 
contribution to the time lag. There is growing evidence that the new branches are formed 
directly at the free barbed ends of the actin filaments, that are in contact with the membrane 
(Carlier,2003). In this case, there would be no other source of time lag. However, if the 
Arp2/3 can diffuse into the bulk and nucleate new branches at barbed ends further back 
from the leading edge, it would be another source of time-lag. Experimental observations 
suggest that barbed ends are localized to a width of A ~ lOOnm from the leading edge 
(Bailly,1999). The average time for a new filament to grow and cover that distance back to 
the plane of the membrane, and add to the pushing force, is: t\ — X/v p ~ 300msec. On this 
time scale, the membrane density fluctuations can diffuse away over length-scales smaller 
than \ft\D ~ 200nm. This length is of the order of 4 unit mesh sizes of the actin gel, while 
we are interested in a continuum description that is valid over longer length scales. We 
therefore conclude that, within these limitations, for length-scales longer than 200nm, we 
can neglect the time-lag. 
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III. RESULTS 

We now solve for the dynamics of the system by first fourier transforming both Eqs. (jll2|) . 
and using solutions of the form: e ~' l ( u ' t +<i- T ) . This gives the following system of equations in 
matrix form 



V 



" 1 | -ifnQ 

h 



(4) 



-iuj + u D —Bq 
—A —iuj + u, 

where B = AkH and ujd — Dq 2 . The dispersion equation of the protein density n and 
membrane height h respectively, is given by equating the determinant of the matrix in Eq0] 
to zero, which yields 

i- (lo d + u q - 4ABq 4 + u 2 D - 2u D u q + ufj 

i~ (uj d + u q + <J AABq 4 + uj 2 d - 2u D u q + ufj (5) 

where, the time-dependent response is given by (n, h)(t) cx exp (—iuj n ^t). It is to be noted 
that these solutions decay exponentially in time if Re[uj nt h) > 0. We now discuss the results 
for two different cases, positive (H > 0) and negative (H < 0) spontaneous curvature of the 
membrane activator proteins. 



—i- 



A. Positive spontaneous curvature 

Membrane activators with positive spontaneous curvature (H > 0) will prefer to aggregate 
at the locations with maximum local curvature (Fig. 2). The solution for the membrane 
height h, in this case, decays with time {Im[ujh\ > 0,Re[uh] = 0), while for the membrane 
density of the activator proteins, n, we find that there can be an instability in the form of 
an exponentially increasing function of time (Im[u n ] < 0,Re[u n ] = 0)(Fig.3a). Depending 
on the form of the membrane response we choose, we get unstable behavior for the following 
range of g-wavevectors 

AriAB 4:T]Vo\H\ 

u q : q < q c , q c = — = 

kD n x 

kD aD 
^ : q > q ^ ^ = ^¥AB ° r 2^45 (6) 

These results arise when the bare response of the protein diffusion is faster than the re- 
sponse of the membrane. In this case, the protein aggregates in response to the mem- 
brane curvature fluctuations, and builds up density fluctuations as it responds faster 
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(Fig. 2). Similar instabilities due to aggregation of membrane proteins with positive 
spontaneous curvature were described in previous studies of different active membranes 
(Chen,2004;Manneville,2001;Ramaswamy,2000). In Fig. 3a we plot u n and Uh as a function 
of q (for the case of a free tension-less membrane). 

The general form of the instability criterion, q < q c (EqEJ), follows from comparing the 
time-scales of membrane motion and in-membrane diffusion of the combined shape-density 
undulations (Fig. 2a). These undulations combine a local increase in the membrane protein 
density, with the driven (active) normal motion of the membrane. The motion of these 
shape-density undulations can be described by an effective diffusion with a dispersion relation 
given by q 2 = Ub U m P /D', where D' = AB/D = Vq\H\k/tiqX- We now discuss the parameters 
that control this motion. The membrane "bump" diffuses faster when the driving velocity 
produced per membrane protein, proportional to Vq/uq, is larger. This is because, density 
fluctuations are converted faster into a height undulation. Larger spontaneous curvature, 
H, causes the membrane bumps to have smaller wavelengths, which move faster. Finally, 
larger osmotic pressure of the membrane proteins \ results in a larger wavelength of the 
density fluctuations, resulting in slower motion. The unstable regime occurs for wavevectors 
where the membrane response, u q , is slower than the rate of diffusion of density-bumps, 
^bump- In this regime, the aggregation of the membrane proteins can occur before membrane 
undulations decay away. The criterion appearing in (jUJ) is simply a restatement of this result. 
It is to be noted that the final expression for q c (EqJHJ) is independent of both the membrane 
bending modulus, k, and the membrane protein diffusion coefficient, D. Quantitatively, 
using the parameters given before (see section II), we find: D' ~ D/2. 

The instability we describe above does not lead to a real divergence, since the velocity 
saturates at v p and the local density of the membrane proteins saturates due to their finite 
size. Note that a similar behavior of the critical wavevector of the membrane instability, was 
shown in (Stephanou,2004): k c oc k a /D a , where k a is the rate of actin polymerization, and 
D a is the bulk diffusion coefficient of actin monomers in the cell cytoplasm. Comparing to 
our expression for q c (EqJHJ), we see that both results are directly proportional to the rate of 
actin growth (k a or vq) and inversely proportional to the diffusion coefficient in the plane of 
the membrane, which tends to smooth away the density accumulation (in (Stephanou,2004) 
the actin diffusion was assumed to be confined to a sub- membrane layer). The details of the 
two models are nevertheless very different. 
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B. Negative spontaneous curvature 

For negative spontaneous curvature of the membrane proteins (H < 0), we find that there 
is a range of wavevectors for which the response frequencies u n and uih are real (Fig. 3b). This 
corresponds to a wave-like behavior, though still damped (or even over-damped) (Fig.2b). 
The range of the wavevectors over which this occurs is given by the condition that the 
square-root term in Eq|3] becomes negative. This gives 



uj q : q < q w , q w = 



K 

The real parts of the frequencies for the height function and the activator density func- 
tion are in anti-phase with each other (Fig. 3b), and correspond to an effective propagation 
velocity v e ff = Re[u]/q in the limit of small wavevectors (q — ► 0) 

1 



u q : v eff = -q^(4D> + D)D (8) 

The damping of these waves is given in this limit by the factor: e~ Dq2t ^ 2 , coming from the 
membrane protein diffusion. For the tension-dominated (second case in Eq|7J), there are 
propagating waves in the limit of large wavevectors (q — > oo), where we get 

1 



<" s> i : v effA = -q^(4D> + D)D (9) 
and the damping of these waves is given in this limit by the factor: e~ ( - Dq2+aq ^ 2v ^ t ^ 2 . The 



waves therefore decay over a length-scale given by: I ~ <7 _1 y (4-D' + D)D/D, which is of the 
order of the wavelength of the density-height perturbation. 

In Fig. 3b we plot the imaginary and real parts of u n and u>h as a function of q (using uj q 
of the flat and tension- less membrane). Note that for q < q w the imaginary parts of both 
frequencies are equal, while the real parts have the same magnitude but opposite signs. For 
any choice of parameters, the motion of n and h changes from damped (i?e[u;] > Jm[o;]) 
to over-damped (i?e[u;] < 7m[o;]) wave-like propagation as q increases (see inset of Fig. 3b). 
Finally, when q > q w , we find the usual exponential decay for both functions (i?e[u;] = 
0,lm[uj} > 0). At the critical wavevector q w , both the response frequencies have the value: 
uj h = u n = -m(VAB + J D)(2 V / AB + Dfifjn 2 . 

In the regime of wave-like propagation q < q w (or q > q Wj i), the density fluctuations 
travel faster than the bare diffusion, due to the additional curvature driving force fEqs l8!9|) 
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(Fig. 3b). The driving force for the wave- like propagation of density-curvature fluctuations is 
shown schematically in Fig. 2b. A local increase in the protein density will result in increased 
membrane curvature there, which then drives these proteins into lower density areas due to 
their negative H, in addition to the usual diffusion. This curvature-induced restoring force 
gives rise to the (albeit damped) oscillatory behavior. 



C. Protein density and membrane height correlations 

So far, we have described the wavevector regimes where, depending on the sign of the 
spontaneous curvature of the activators, one gets either an instability or wave-like modes. 
We now wish to address the question of the actual amplitudes of the fluctuations that 
characterize the motion in these regimes. This we do by calculating the correlation functions 
of the membrane height and of the membrane activator density. Solving Eq0]we get 

, a , ,v (f 2 (g,u;))q 2 (uj 2 + uj 2 ) 

{Tl W ' U)) u\u D + uu q y + (u D u q - u; 2 - ABqtf 

, h 2( U (/n(g»^))g 2 ^ 2 / 1n N 



Integrating these functions over uj we find the spatial (static) correlations (n 2 (q)), (h 2 (q)). 
We plot these functions in Fig. 4, for a free and tension-less membrane. Analytic expressions 
for these functions can be calculated, but are quite lengthy, so we will give them explicitly 
only for the limiting cases. 

It is to be noted that the height fluctuations in our model are derived solely from 
thermally-driven density fluctuations of the membrane proteins. This is why we get: 
{h 2 (q,u)) oc (f 2 (q,ou)) oc k B T. These height undulations are superimposed over the av- 
erage forward motion of the membrane, at average velocity vq, determined by the average 
density n . 

Another notable point is that we assumed a continuous and constant force, or driving 
velocity, due to the actin polymerization in Eq|2] This is reasonable as long as we are 
interested in membrane motions on timescales longer than the duration of an individual 
actin polymerization event. More generally, we can describe the actin-induced velocity (or 
force) by a random shot-noise behavior (Gov,2004), with a typical time r. This amounts 
to replacing: A 2 —>■ A 2 /(l + (ujt) 2 ) in the numerator of Eq^Jfor (h 2 (q,u))). This is easily 
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calculable, and is found to change the behavior quantitatively, but not to change the value 
of the critical wavevector q c , or the qualitative forms of (h 2 (q)) in the q — > 0, oo limits. 

We now discuss the form of the density and height correlations in various cases and limits. 

For positive spontaneous curvature of the membrane proteins (H > 0) (Fig. 4a), we find a 
divergence of both the density and membrane height fluctuations at the critical wavevector 
Qc (Eq. Around the critical wavevector, setting q = q c + 5, the divergences have the form 



2 , _ (f 2 ) <£« _ IrrvlH k B Tn 1 



(n\q c )) 



2D At]\5\(D + D') x 3 M D + D ') l<*| 



(h 2 (Q ^ = MPJA 2 = k B T irrvj 1 

^ {qc,/ V2q^\5\D(D + D') Kqi V2 X n (D + D')\5\ 1 } 

where, from the last line, we can define: (h 2 (q c )) = ksT e ff/Kq^. The "effective- 
temperature": T e ff/T = Ar]v 2 /\/2xn (D + D')S, diverges at the critical wavevector. This is 
reminiscent of the divergence in the effective temperature describing the response of active 
hair bundles in the hair-cells of the auditory system (Martin,2001) when there is a resonance 
with an internal driving frequency. In our case the divergence occurs when the length-scale 
of the active process is in "resonance" with the length-scale of the spontaneous curvature. 

In the limit q — > 0, the height correlations have the form: (h 2 (q)) = kBT e ff/K,q A . Here 
we chose to define an effective temperature, T e ff, since the power law dependence is similar 
to the behavior of the thermal membrane height fluctuations (Safran,1994): (h 2 (q))T = 
ksT/ nq A . The appearance of thermal-like correlations is not surprising, since the driving 
force for the height fluctuations comes from the thermal fluctuations of the membrane protein 
density: {f 2 (q, u>)) oc ksT. Thermal-like correlations also appear for various choices of active 
membranes (Manneville,20001;Ramaswamy,2000;Gov,2004). The effective temperature we 
defined, has the following limiting forms 



D' -> 



T eff v KV 2 V 



T 2DD'n oX 2D\H\ 



D'^oo ^ - KV ° = - ! V ° n ^ (12) 
T 2v / ^D 7 % oX 2y Dk\H\ 3 V ; 

The functional form of T e ff/T is very intuitive: the effective temperature increases with the 
pushing velocity of the actin vq, and is inversely proportional to the diffusion coefficient of 
the membrane proteins, which smooths away the density fluctuations. 

In contrast, the density fluctuations are finite in the limit q —* 0. Note that for the free 
diffusion problem, we recover the usual free diffusion: (n 2 (q)) = n q 2 k B TD/2w D ~ n . In 
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our model we find in the following limits 

D' -> <n 2 (0)> -> I i 

X 

D^oo (n (0)) -» - y— (13) 

The first limit shows the approach to the bare membrane diffusion in the absence of actin 
polymerization (vq — > 0). In the second limit we find that the rapid formation of membrane 
undulations, due to t; oo, effectively localizes the membrane proteins and suppresses any 
long wavelength density fluctuations. 

In the limit q —>■ oo the density and height fluctuations are given by 

2/ u (fn) n k B T 



(n 2 (oo)) 



2D X 



= 16^1 (14) 
X n o K Q 

The density fluctuations are finite, and approach the bare membrane diffusion result (see 
first part of EqfT3*j). The height fluctuations decay in this limit in a non-thermal form, 
reminiscent of similar results for model active membranes (Gov,2004). 

In the case of negative spontaneous curvature of the membrane proteins (H < 0), as 
we have already seen, no instability occurs (fig.4b). For the density correlations we find 
that (n 2 (q)) is approximately a constant as a function of q (Fig.4b), close to the limiting 
value (n 2 (oo)) (Eqll4|) of free diffusion. There is a region of reduced correlations, that 
corresponds to the additional restoring curvature force, which now acts to smooth away any 
density fluctuations. The density correlations therefore dip for wavevectors q < q w (EqJTJ), 
where there is wave-like propagation. The height correlations show a cross-over from 1/q 4 
to 1/q 6 decay around q w . 

We now consider the correlations for the tension dominated case (or thin lamellipodium 
edge). Here the membrane response is given by oo qt i,u qjt oc q. The results in the low and 
high q limits are 

2 k B T eff T eff Attv^t] 2 

q^O (h 2 {q)) ~ f- = — 

aq z 1 n x<J 

n 2 I \\ ^BT e ff 

q^oo (h\q)) ~ f- 

nq* 

T VQ ^{2D'/D + 1) + y/W/D + 1 - sj(2iy/D + 1) - ^AD'/D + 1 



T 2D\H\ 



AD' /D + 1 
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(15) 

The first limit has the same form as for thermal fluctuations in the tension dominated 
regime, while the large q limit has the form of thermal fluctuations in a free and tension-less 
membrane, prompting us to express both in terms of an effective temperature. 

One of the most striking results is that the height fluctuations increase with increasing 
fluid viscosity r\ (Eqs JlHfl4pi5|) .This behavior can be intuitively understood as arising from 
a height fluctuation that is driven at a constant velocity vq over a time interval that is 
proportional to the viscosity (l/w q oc 77). This is similar to the results of a previous analysis 
of active membranes driven by ion pumps (Prost,1998). In both these cases, the active 
membrane proteins are allowed to have density fluctuations, while they impose a given 
velocity on the membrane. This is in contrast to the results of active membrane proteins 
that produce a fluctuating force with zero average value (shot-noise), where the height 
fluctuations are found to decrease with increasing fluid viscosity (Gov,2004). Note that if 
the actin imposes a force, F actin , rather than a velocity v , on the membrane, this introduces 
another factor of l/q 2 rf to the height correlations (h 2 (q)) (see discussion following Eq|2J), 
and the dependence on the viscosity is subsequently modified. 

The density and height fluctuations also increase when the membrane diffusion coefficient 
is decreased (Eqs ll 1I12I15|) : (h 2 (q)) oc 1/D. A similar result appeared in (Prost,1998), but 
was limited there to the tension-less regime. We find that this behavior also appears in the 
tension-dominated case (EqJ15p. which seems to be more realistic for living cells (see next 
section). The origin of this behavior is very intuitive; slower diffusion allows fluctuations in 
the membrane protein density to survive longer, causing larger height undulations. 

We now examine how the correlations depend on frequency. The power spectrum of the 
height fluctuations as a function of frequency can be obtained by integrating Eq^Jover q. 
This gives us the temporal correlations (n 2 (uj)), (h 2 (uj)). Since the integration is not possible 
analytically, we do it numerically. We plot the height correlation function in Fig.5, for the 
cell membrane with the elastic parameters: k and a. These parameters are found empirically, 
by fitting to the observed thermal fluctuations alone, i.e. when the actin polymerization is 
blocked (Zidovska,2003) (see next section). For comparison we also plot the correlation that 
arises purely from thermal fluctuations. This approaches the free membrane limit at high 
frequencies u — > 00, where (h 2 (u)} thermal —> uj~ 5 ^ 3 (Zilman,2002). At lower frequencies, in 
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the tension dominated regime, the thermal behavior is: oc uj^ 1 . 

The active fluctuations are found to have a {h 2 (u)) act i n oc u~ 2 behavior at small uj — > 
0, and (h 2 (uj)) actin oc uj' 3 at large uj — > oo. The cross-over occurs roughly where the 
frequency of the membrane bending modes equals the frequency of the effective diffusion 
of the membrane bumps. For the case when H < 0, the cross-over corresponds to the 
appearance of propagating waves and occurs at the wavevector, q w . The frequency of the 
height fluctuations corresponding to this wavevector, is given by u>h, shown as the vertical 
dashed-line in Fig. 5. 

Future experiments that probe the power spectrum of the height fluctuations, should show 
a clear difference between the thermal and active components. The actin-induced contribu- 
tion to the height fluctuations are predicted to be much more confined to low frequencies 
than the thermal contribution. 

IV. DISCUSSION 

A. Comparison with experiments 

Experimental observations of the actin-driven motion of cells and cell membranes 
(Gerisch, 2004; Vicker, 2002; Giannone, 2004) show both wave-like propagations and finger-like 
filopodia. Since our model predicts that both behaviors are possible, depending on the sign 
of the spontaneous curvature, H, of the membrane proteins, one possibility is that activators 
with both types of spontaneous curvature exist in vivo. There is also the possibility that 
the spontaneous curvature of a membrane activator is altered by a conformational change 
that is brought about by phosphorylation or binding to another protein (or a number of pro- 
teins), either already in the membrane or from the cytoplasm (Small, 2002). Thus the cell 
has many options, all of which it may use to adjust the local concentration of the membrane 
proteins which produce either uniform growth or filopodia (Biyasheva,2004;Nozumi,2003). 
We predict that these proteins (or protein complexes) have different spontaneous curva- 
tures: H < for uniform growth and H > for filopodia. Indeed for filopodia growth 
the membrane proteins have to form specific complexes which, in light of our model, must 
have large spontaneous curvature (Gauthier-Campbell, 2004; Wood, 2002). These complexes 
can then recruit crosslinking proteins such as fascin, which help form tight actin bundles 
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inside the growing filopodia (Biyasheva,2004). Future extensions of our model may include 
the coupling of dynamical changes in the spontaneous curvature, H, to the local densities 
of various proteins. 

Another source of experimental corroboration for our overall qualitative picture of the 
dynamic co-localization of activator membrane proteins and actin polymerization in high (or 
low) curvature regions of the leading edge, is provided by the experiments of Nalbant et al 
(Nalbant,2004). The authors visualized the dynamics of activated Cdc42 in living cells and 
showed that there was a strong correlation between the most recently formed protrusions and 
the level of active Cdc42 in its vicinity, mostly concentrated at the tips of the protrusions. 
They also show that the activator proteins are present only near the base regions of filopodia, 
consistent with our assertion that the activators form the initial seed for the filopodium by 
assembling an actin rich bump, whereafter the newly recruited bundling agents and normal 
barbed end polymerization of the actin filaments can lead to the filopodia structure. 

For positive spontaneous curvature of the membrane proteins (H > 0), we can estimate 
the critical wavevectors (EqEJ), using typical values of the various parameters (see sec. II). For 
the spontaneous curvature we use H ~ (5 — 100nm) _1 (Girard,2005). For the flat membrane 
case we find: q^ 1 ~ 1 — 10/im. For the membrane edge case we find: q~l ~ d ■ (dq c ), 
where typically: d ~ 0.5 — 1/im. Our analysis predicts a specific wave-vector (q c ) which 
becomes unstable, so that the resulting filopodia should have an average spacing given by 
the corresponding wavelength. This length-scale appears to correlate well with the observed 
average separation between neighboring filopodia (Oldenbourg,2000), of 1.5 — 3/im. Note 
that from EqlHl increase in the membrane tension causes an increase in the density of 
filopodia (Parker, 2002). Increased membrane tension was found to reduce the velocity with 
which actin polymerization is pushing the membrane (Raucher,2000), so that we expect not 
only more numerous but also smaller filopodia under increased membrane tension. 

We predict that the density and height fluctuations increase when the membrane dif- 
fusion coefficient is decreased (Eqs Jllll2ll5j) . The membrane diffusion coefficient may be 
changed by addition of various chemical agents, such as changing the cholesterol content 
( Vasanji, 2004; Vrljic, 2005). Note that changing the cholesterol level may affect the activa- 
tion of the membrane proteins (Niggli,2004), which is a process that is not included in the 
present work. 

This prediction may explain the recently observed low membrane diffusion (high micro- 
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viscosity) at the leading edge of moving cells (Vasanji,2004). The ratio of the diffusion 
coefficients between the cell side and leading edge is found to be: D trai i/ Di ead ~ 3. If we 
correlate the r.m.s. height fluctuations from Eqs. (jllJ12j) with the average rate of lamel- 
lipodial extension, we predict: Vi ead /V side ~ sj (h 2 (q)) lead / (h 2 (q)) side ~ ^ D trail / D lead ~ 
This is in good agreement with the measurement (Vasanji,2004). Our model therefore pro- 
vides a natural explanation for this otherwise paradoxical observation: the membrane is 
stiffer (more viscous) in regions where motility is increased. Presumably larger undulations 
in the shape of the leading edge, help the cell overcome local friction barriers, and results in 
faster overall motion (Ehrengruber,1996). Additionally, the membrane undulations at the 
moving front can provide localization points for the formation of adhesion complexes, which 
are important in completing the cycle of cell motility (Lavelin, 2005; Zaidel-Bar, 2003). 

This result of our model may also explain the observed response of endothelial cells' 
motion to shear flow (Tardy, 1997; Albuquerque, 2003). In these experiments it was shown 
that the cells move less quickly against the direction of the flow, as compared to the per- 
pendicular and parallel directions. Concurrent with this motion, there is an increase in 
the fluidity of the membrane in the front part of the cell, by as much as a factor of 2 
(Haidekker,2000;Butler,2001). According to our model the amplitude of the active mem- 
brane fluctuations is therefore reduced by a factor of ~ v2 compared to the rear of the 
cell, which is in very good agreement with the measured drop in the fluctuation amplitude 
in the presence of shear (Dieterich,2000). This then results in the observed orientational 
motility (Tardy, 1997; Albuquerque, 2003). Such a physical response to shear flow may com- 
plement or trigger the biochemical changes that take place in the presence of shear-flow 
(Zaidel-Bar, 2005). 

For negative spontaneous curvature of the membrane proteins (H < 0), we can esti- 
mate the critical wavevector (EqJJJ) using the values of the parameters that appear above. 
This gives us q^ 1 ~ 1 — 10/im, which is similar to what we obtained for q c (EqlHJ). Thus 
both the instabilities and the wave-like motions have the property that they occur only 
for membrane length-scales larger than some critical length-scale, ~ 1/zm. Indeed, there 
are no long-lived actin structures smaller than this length-scale, on the cell membrane 
(Gerisch,2004;Vicker,2002). 

For the velocity of the propagating waves, we get from EqlHJ v e ff ~ qD, which results 
in velocities of the order: v e ff ~ 1 — 0.1/xm/sec for wavevectors: q ~ l(fim)^ 1 . Actin 
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waves with these velocities and wavelengths are indeed observed on the surface of cells 
and lamellipodia (Gerisch, 2004; Vicker,2002;Giannone, 2004). Note that in these experiments 
the observed waves are on the bottom part of the cell, where the membrane is largely 
flat next to the glass substrate. Our analysis predicts that the actin waves correspond to 
small undulations on the membrane surface (Fig.2b). Our interpretation of these waves is 
therefore different from that given in (Vicker,2002), where the surface waves are proposed 
to be sections through three-dimensional spiral waves in the cell bulk. Recent experiments 
seem to confirm our interpretation since they suggest that the actin structures are largely 
confined to the cell membrane (Gerisch, 2004), and that the travelling- wave propagation is 
related to actin polymerization around a high density Arp2/3 complex (Bretschneider,2004). 
The formation and decay of these density fluctuations occurs on a timescale, tfi uct ~ 2 — 
3sec, which is in agreement with membrane diffusion times over the lateral size of these 
formations (~ 1/2/im). The slower decay as compared to the formation, ~ 3.5sec vs. 
~ 2sec respectively, may be due to the extra distance to diffuse out of the bump (Fig.2b): 
tfiuct ■ vo ~ O.lyum. Alternatively the diffusion coefficient may decrease due to the dense 
actin gel formation. 

Recently the static height correlations, (h 2 (q)), were measured on living cells (Zi- 
dovska,2003). The mean-square height undulations of the active cell are found to be ~ 8 
times larger than for the inactive cell (Zidovska,2003). The correlations agree with the 
tension-dominated behavior given in EqE3 (Fig.6), if we use the same parameters as before 
and take the surface tension to be a ~ 0.5 x 10~ 8 J/m 2 . In particular, cells that lack the 
actin-polymerization motility, display much smaller fluctuations (Zidovska,2003) presumably 
of thermal origin. Indeed these fluctuations are well described by confined thermal correla- 
tions of the form (Gov,2003): (h 2 (q)) oc k B T/{Kq A + oq 2 + 7), with a ~ 3.7 x 10- 8 J/m 2 and 
7 ~ 2.6 x 10 5 J/m 4 (Fig.6). The equivalent confinement distance: (It — k B T /R^ffR, ~ 70nm, 
is consistent with the average separation between the fibers of the actin mesh, which under- 
lies the membrane. In a normal cell, when actin polymerization is driving the membrane 
fluctuations, there is no meaning to any membrane confinement. It is then possible to fit the 
active cell data to the thermal fluctuations of an unconfined membrane, with an "effective 
temperature" of T e ff/T ~ 8 (Fig.6). This approach though does not give us any informa- 
tion about the nature of the active fluctuations. From the q — > limit of EqJTol we get an 
effective temperature of this magnitude if we take: a to be smaller than the value given by 
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the fit to the thermal fluctuations, by a factor of ~ 8, n ~ (300nm) -2 , and an effective 
viscosity 77 ~ 100r] water - These parameters are within reasonable limits for a cell, but an 
exact comparison with the data awaits independent determination of these parameters. 

Other experimental data (Agero,2003) indicate that in a living cell there is an exponential 
component to the probability distribution function of the height fluctuations, while in a cell 
that has its actin polymerization blocked, only the usual Gaussian term due to thermal 
membrane fluctuations remains. We can assume that the usual thermal height fluctuations 
add incoherently to the active fluctuations, which according to our model arise from the 
fluctuations in the density of the membrane proteins, n. For simplicity let us describe 
the case where H = 0, so that the density of the membrane proteins is decoupled from 
the membrane height fluctuations (Eq^). In the limit of small density fluctuations, the 
corresponding actin polymerization induced height fluctuations are given by EqJ2 

AT? 

5h ~ A5ndt 5n 2 (16) 

where the diffusion occurs in a membrane patch of size L over a time: dt ~ (L 2 /D)(5n/no). 
This relation states that larger density fluctuations (larger 5n/no) take longer to diffuse 
away. The probability distribution function of the density fluctuations Sn/n Q is a thermal 
Gaussian (Manneville,2001;Ramaswamy,2000), so that combined with EqJ16l we get 

P(5n) ~ e (-xSnVnlk B T) 

P(5h) ~ e (-xDSh/v () L^k B T) ( 17 ) 

Thus the probability distribution does indeed have an exponential component. The larger 
the driving velocity Vq, the more enhanced the distribution for larger height fluctuations. On 
the other hand large values of the diffusion coefficient D or membrane protein compression 
energy x, narrow the distribution by causing density fluctuations to have shorter life-times 
and amplitudes. A quantitative comparison with the observations awaits future experiments. 

Recent data (Neto,2005) indicates that the velocity of actin induced membrane ruffles, 
is strongly temperature dependent. According to our model this velocity is proportional to 
the membrane diffusion coefficient fEq J8!9|) : v e ff oc Dq, where q is the inverse of the typical 
length-scale of these ruffles (usually ~ 1 — 2/im (Neto,2005)). Since the diffusion coefficient 
is inversely proportional to the membrane viscosity (Almeida, 1995), we expect it to vanish 
at the liquid-gel transition temperature T m , where the viscosity diverges (Dimova,2000): 
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D(T) oc \T— T m \ 1A . Using typical values for T m = 20°C, and q = 0.67 ^m -1 , we fit the overall 
scale of D(T) (inset of Fig. 7), such that the resulting velocity v e ff agrees with the observation 
(Fig. 7a). Furthermore, the observed decay time of smaller height fluctuations, also behaves 
as: Tdecay oc ^-/Dq 2 , with a larger wavevector (smaller wavelength) of q = 1.15/im -1 (Fig. 7b). 

Finally, the mean-square membrane curvature was observed not to depend on the temper- 
ature (Neto,2005). The mean-square curvature: (H 2 ) = J q 4 h 2 d 2 q, is dominated by q — > oo 
modes, while the mean-square amplitude of height fluctuations (h 2 ) is dominated by the 
q — > modes. From our model we predict that for the tensionless membrane, the amplitude 
of the q — > modes does depend on the diffusion D (EqEJ): (h 2 ) oc l/D(T), while the 
amplitude of the q — > oo modes does not (EqJ14jl. This difference could explain the inde- 
pendence of the observed r.m.s. curvature on the temperature, except for the very weak 



(H 2 ) oc \fksT which amounts to ~ 2% over the observed temperature range (Neto,2005). 



Detailed comparisons between our model (namely (h 2 (q)) ,(h 2 (uj))) and the experimental 
data (Gerisch,2004;Bretschneider,2004), awaits more quantitative analysis of the spatial and 
temporal shape fluctuations in living cells (Zidovska,2003;Agero,2003;Neto,2005). 

So far we have discussed the membrane dynamics at the leading edge which is what we 
are modeling. However, our model can also give us insights into phenomena that occur at 
the cellular scale. An example of complicated, oscillatory dynamics of the bulk actin gel, 
is described in (Giannone,2004). The authors find a periodic interruption (~ 20s) in the 
forward motion of lamellipodia that had no filopodia. One possibility for the mechanism is 
shown schematically in (Fig.8). The activating membrane proteins are initially concentrated 
at the lamellipodium edge, and since there are no filopodia we can take the spontaneous 
curvature to be small or negative. As the forward motion persists the leading edge thins 
such that the local curvature is too high for these proteins and they prefer to move to the less 
curved membrane on the upper surface. This causes a backward propagating wave of actin 
polymerization, which proceeds until the back edge of the lamellipodium. This explains why 
the contractions occur every t cont ~ Li am /v ~ 10 — 30sec, where Li am ~ 2[im is the thickness 
of the lamellipodium (Giannone,2004) and v ~ 0.1/xm/sec is of the order of the calculated 
propagation velocity v e ff (Section IV) (Fig.8). The membrane dynamics we considered in 
the present model are therefore coupled in the cell to the dynamical variations affecting the 
entire actin layer, and this coupling remains to be described. 
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B. Predictions 

Our model allows us to make testable and quantitative predictions. For example: 

• Changing the fluid viscosity and membrane tension will shift the average density of 
filopodia (EqEJ). 

• The velocity of propagation of actin density fluctuations on the cell membrane, is 
predicted to be linear in the diffusion coefficient of membrane proteins (EqJSJ). 

• Similarly, the amplitude of density and height fluctuations are predicted to increase 
when the membrane diffusion coefficient is decreased (Vasanji,2004) (Eqs Jllll2ll5jl . 
Note that the density of filopodia, given by q c (JUJ), is independent on D. 

• In the tension-dominated regime, which seems to be applicable to most cells, the long 
wave-length height fluctuations increase with increasing of the fluid viscosity (EqJISp. 
The same behavior is found also in the tension- less case (Eqll2|). 

• Our prediction that the membrane proteins that initiate filopodia (Nozumi,2003) have 
a high spontaneous curvature, has to be tested. Incorporation of these proteins into 
synthetic vesicles and observing the resulting shape transformations could determine 
this parameter. 

Some of these manipulations are possible in living cells, while others are better tested in 
synthetic systems (Vasanji,2004). 

V. CONCLUSION 

The dramatic membrane dynamics that occur at the surface of stimulated cells is a 
consequence, not only of the actin polymerization dynamics, but also of the interplay between 
the dynamics of the membrane itself and that of the activators that reside on it. Keeping 
this in mind, we presented a simple model that treats the dynamics of a membrane under 
the action of actin polymerization forces that depend on the local density of freely diffusing 
activators on the membrane. We took into account the thermal density fluctuations and the 
spontaneous curvature associated with the activators and showed that, depending on the 
spontaneous membrane curvature associated with the activators, the resulting membrane 
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motion can be wave-like, corresponding to membrane ruffling and actin-waves, or unstable, 
indicating the tendency of filopodia to form. Thus our simple model system managed to 
capture the wide range of complex dynamics observed at the leading edges of motile cells 
both qualitatively and quantitatively indicating that the essential physics had been retained. 
Our model not only provides detailed estimates of the morphology and dynamics of the 
membrane structures, but also provides quantitative explanations for a variety of related 
experimental observations. These include the puzzling increase in membrane micro- viscosity 
at the leading edge of migrating cells, the appearance of an exponential contribution to the 
probability distribution of membrane height fluctuations and the temperature dependence 
of the membrane ruffle velocity among others. Thus, our model offers a simple framework 
with which to analyze and understand experimental data and make quantitative predictions 
to be tested by future experiments. 

We should, however, keep in mind that cell motility involves many processes which we 
did not take into account in our model, such as adhesion, formation of stress fibers and the 
action of molecular motors. Even within the context of our model, the dynamics of the actin- 
driven cell motility is largely assumed to be controlled by dynamics of proteins on the cell 
membrane, which fails to capture the link between the dynamics of these membrane proteins 
and the bulk dynamics of the actin gel, occurring behind the moving front (Plastino,2004). 
Integrating all these components into a holistic picture remains a challenge. We should 
therefore view our model as representing the physical dynamics of the membrane-actin 
system, which trigger the formation of patterns in the membrane morphology, that are 
part of the overall motility mechanism. Our model therefore provides answers to one part 
of the overall problem of cell motility and should be useful for any integrated approach to 
cellular motility. 
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VII. FIGURE CAPTIONS 

Fig.l: Schematic picture of the model, (a) The Arp activating membrane proteins are 
symbolized by the black squares, diffusing in the flat membrane. Where they have a high 
density the actin polymerization is more extensive (dashed regions) and so is the velocity of 
the membrane (normal arrows), (b) In the case where the polymerization is confined to a 
thin leading edge, the local high curvature changes the response to co> g l . 

Fig.2: Schematic picture of the two behaviors depending on the spontaneous curvature 
of the membrane proteins, (a) H > 0: fluctuations in the density of the Arp activat- 
ing membrane proteins grow unstable when the proteins aggregate into the high curvature 
"filopodia". (b) H < 0: wave-like propagation due to the restoring force of the curvature, 
breaking up high density fluctuations (dashed arrows). 

Fig. 3: Calculated response frequencies of the membrane protein density u n and membrane 
height undulations Uh (using oj q of a flat membrane, k = 10k B T and v = 1/xm/sec). The 
bare diffusion rate ujj and bare membrane response uj q are shown by the dotted and dot-dash 
line respectively, (a) H = (10nm) -1 : co^-solid line, c^-dashed line. The critical wavevector 
q c (EqJOJ) below which the density fluctuations are unstable is indicated by the vertical dotted 
line, (b) H = — (3nm) _1 : The imaginary parts are given by the solid lines while the real 
parts are given by the dashed lines. The critical wavevector q w (Eq|7j) below which the 
wave-like fluctuations occur is indicated by the vertical dotted line. In the inset we show 
that the imaginary part can be smaller than the real part for small enough q. 

Fig. 4: Calculated static density and membrane height correlation functions: (n 2 (g))-solid 
line,(/i 2 (g))-dashed line, (a) H = (lOnrn) -1 : The correlations diverge at the critical wavevec- 
tor q c (EqJHl), indicated by the vertical dotted line. The limiting values (n 2 (0)) (Eq J13Jl and 
(n 2 (oo)) (EqJ14|) are shown by the horizontal dotted and dash-dot line respectively, (b) 
H = — (3nm) _1 : The density correlations dip for wavevectors q < q w (Eq[7J), indicated by 
the vertical dotted line. The height correlations show a monotonous decay, having a cross- 
over from l/q 4 to 1/g 6 behavior around q w . The limiting value (n 2 (oo)) (Eqll4j) is shown 
by the horizontal dash-dot line. 

Fig.5: Calculated height fluctuations as a function of the frequency u. The dash-dot line 
is the thermal fluctuations, that approach lu~ 5 / 3 at large frequencies (asymptotic dashed- 
line). The actin-induced fluctuations are given by the solid and dotted lines, for H > and 
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H < respectively. The asymptotic behavior is given by the dashed straight lines: lu~ 2 
and u>~ 3 in the limit of small and large frequencies respectively The vertical dashed line 
represents the frequency of the cross-over, roughly given by u)h(q w ). 

Fig. 6: Calculated static height correlation function (h 2 (q)) in the tension-dominated 
regime (Eqll5j) (dashed line), compared with the data (Zidovska,2003) for normal cell (stars) 
and inactivated cell (squares). The solid line and dash-dot lines show the behavior for the 
thermal height fluctuations in an unconfined and confined membrane respectively. The lower 
panel shows the two cases: (a) Inactivated cell with confined thermal membrane fluctuations, 
and (b) actin-induced fluctuations in the normal cell. 

Fig. 7: (a) Calculated velocity of actin-induced membrane ruffles: v e ff oc D(T) (solid line), 
where D{T) is shown in the inset. The experimental data (circles) is from (Neto,2005). (b) 
Calculated decay time of actin-induced membrane fluctuations: r decay oc 1/D(T) (solid line), 
compared to the experimental data (squares) (Neto,2005). 

Fig.8: Schematic description of the periodic contractions found in lamellipodia growth 
(Giannone,2004). (a) The membrane proteins initially localized at the lamellipodia's edge 
(filled circles) produce a forward pushing actin network (shaded ellipse). (b,c) As the edge 
thins they are pushed towards less curved regions at the top of the lamellipodia, and then 
propagate until the back edge, (d) The forward edge is again thicker now, so these membrane 
proteins can localize there and forward motion resumes. 
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